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We compute the electronic structure, momentum resolved spectral function and optical conduc- 
tivity of the new superconductor LaOi-ccF^FeAs within the combination of the Density functional 
theory and the Dynamical Mean Field Theory. We find that the compound in the normal state 
is a strongly correlated metal and the parent compound is a bad metal at the verge of the metal 
insulator transition. We argue that the superconductivity is not phonon mediated. 

PACS numbers: 71.27.+a,71.30.+h 



In the Bardeen-Cooper-Schrieffer theory of supercon- 
ductivity, electrons form Cooper pairs through an inter- 
action mediated by vibrations of the crystal. Like lat- 
tice vibrations, antiferromagnetic fluctuations can also 
produce an attractive interaction creating Cooper pairs, 
though with spin and angular momentum properties dif- 
ferent from those of conventional superconductors. Such 
interactions was implicated for class of heavy fermion 
materials based on Ce, U and Pu with rather low transi- 
tion temperatures, and cuprate superconductors with the 
highest known transition temperatures. Recently a sur- 
prising discovery of superconductivity in iron-based com- 
pound LaOi-^F^FeP [1 with T c ~ 7K sparked a new di- 
rection to explore superconductivity in a completely new 
class of materials. Very recently a substitution of P by 
As raised T c to 26 K [2] becoming already one of the su- 
perconductors with highest T c among non-cuprate based 
materials. Exploring superconductivity in similar iron 
based compounds holds a lot of promise for increasing 
T 

In this letter we explore the electronic structure 
and optical properties of LaOi-^F^FeAs within Density 
Functional Theory (DFT) and Dynamical Mean Field 
Theory (DMFT). 

LaOFeAs has a layered tetragonal crystal structure 
shown in Fig|l] Layers of La and O are sandwiched 
between layers of Fe and As. The Fe atoms form a 
square two dimensional lattice with Fe-Fe lattice spacing 
of 2.853A. To understand the material properties, it 
is important to identify the character of dominant bands 
near the Fermi level, their energy and momentum dis- 
tribution. For this purpose, the first principles density 
functional theory is the invaluable tool. We used the 
full-potential augmented plane-wave method as imple- 
mented in Wien code [3 . The lattice parameters and in- 
ternal atomic positions have been determine experimen- 
tally (a = 4.035, c = 8.741, z La = 0.142, z As = 0.651). 
For the exchange correlation potential we used gradient 
approximation [4 (GGA) in the Perdew-Burke-Ernzerhof 
variant, and 12 x 12 x 5 k-points. 

This method predicts that the dominant states at the 
Fermi level come from Fe 3d atomic states extending 
roughly between -2 eV and 2 eV as shown in Fig. [2j 




FIG. 1: Crystal structure of LaOFeAs. The middle layer 
consists of Fe atoms (red spheres) which form a two dimen- 
sional square lattice and As atoms (yellow) . The second plane 
consist of O (green) and La (gray) atoms. To describe the or- 
bital character in Fig. [3] the coordinate system is chosen such 
that both x and y axis point from each Fe atom towards the 
nearest neighbor Fe atoms. 
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FIG. 2: Density of states within the GGA approximation. 
The partial character of Fe, O and As is shown separately. 
Due to presence of As, a lot of electronic charge is found in 
the interstitial regions and can not be assigned an atomic 
character. 



The important feature of LaOFeAs compound is that 
DFT predicts a very steep and negative slope of the den- 
sity of states (DOS) at the Fermi level. In the rigid 
band approximation, neglecting the many body effects, 
the hole doping will lead to increases of the DOS at the 
Fermi level, while the electron doping will be accompa- 
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nied by a decrease of DOS. Hence, for the conventional 
phonon mediated superconductor, one would expect de- 
crease (increase) of T c when doping electrons (holes), con- 
trary to what is observed in experiment [2]. According 
to Ref. [2], the electron doping by F~ leads to super- 
conducting ground state, while Ca 2+ hole doping does 
not cause superconductivity. This gives a clear hint that 
the superconductivity in this compound might not be 
phonon mediated. Indeed an explicit calculation of the 
phonon coupling constants within the DFT, gives a value 
too small to explain the observed critical temperature. 

While the dominant electronic character near the 
Fermi level is due to Fe, a strong mixing with As is ap- 
parent at — 2.7eV where As 4p band is strongly peaked. 
There is a noticeable similarity between the Fe-As mixing 
and oxygen-transition metal mixing in transition metal 
oxides. However, there is an important difference be- 
tween the two. The As bands are much more itinerant 
and broader than oxygen 2p bands. Consequently, plenty 
of As electronic charge is found in the interstitial region 
far from any of the constituting atoms, which is shown 
separately in Fig. [2] 

While the crystal field splitting in most of transition 
metal oxides leads to a clear separation of transition 
metal d bands into ti g and e g part, this is not the case 
in LaOFeAs due to the itinerant character of As. All five 
Fe-d bands therefore participate in the bonding and in 
the Fermi surface of the compound. 

Within DFT, the parent compound is very close to 
magnetic instability and at normal pressure DFT selects 
ferromagnetic state while 10% compressed solid preferes 
paramagnetic state. The magnetic moment predicted by 
DFT is very small of the order of 0.15/i^. We note in 
passing that within DFT the antiferromagnetic state is 
unstable at normal pressure. 

The bandwidth of iron 3d bands is only around 3eV 
while the Coulomb repulsion in this orbital is around 4 eV 
[5]. For a single band system with this bandwidth and 
Coulomb repulsion, one would expect a Mott insulating 
ground state. On the other hand, the perfectly degen- 
erate system with five degenerate d bands would still be 
below the Mott transition, since the critical U for the 
disappearance of the Fermi liquid solution within DMFT 
is proportional to AT, where N is the band degeneracy of 
the system (in case of d electrons 5). The five d bands 
of Fe in LaOFeAs are not perfectly degenerate, but one 
would still expect relatively large critical U for the Mott 
transition for the odd number of electrons in the d band. 
In the parent compound LaOFeAs, the number of d elec- 
trons is 6. The Coulomb repulsion can strongly enhance 
the splitting between the bands [6] and the gap can open 
at the Fermi level. This interplay of crystal field split- 
ting and correlation effects was addressed in many model 
Hamiltonian studies [6j [TJ [8] and it was shown how a bad 
metal or bad semiconductor can appear on the metallic 
side at finite temperature. In LaOFeAs we checked that a 



slightly enhanced Coulomb repulsion (U = 4.5 eV) leads 
to a finite gap in the d band. For typical Coulomb re- 
pulsion of Fe (U = 4eV) the system is still metallic, 
but a bad metal having, some characteristics of a bad 
semiconductor. The metallic state is however still very 
correlated with quasiparticle renormalization amplitude 
between Z ~ 0.2 - 0.3. 

To describe this type of system, one needs to go beyond 
the traditional band structure methods and concentrate 
on the spectral function A(k, uj) ll> of the system which 
described the probability to add or remove an electron 
with momentum k, angular momentum L = (/,m) and 
energy uj. 

To compute the spectral function of the LaOFeAs, 
we used Dynamical Mean Field Theory (DMFT) [10] 
which takes into account the strong Coulomb repulsion 
among the correlated set of bands and can describe the 
dual nature of the electrons in correlated compounds; the 
itinerant part of spectra which forms narrow quasiparti- 
cle bands at low energy, and the localized part of spectra 
at higher energy, which correspond to the nearly localized 
atomic nature of the electron. In this method, the spec- 
tral function is obtained from the one electron Green's 
function A(k, uo) = (G^ (k, uj) — G(k, uj))/ (2iri) where the 
latter takes the form 

G(k ' w) = o k (. + M )-F k -EH- (1) 

The one electron part of the Hamiltonian and overlap 
matrix Ok is obtained by the LDA method [11 while the 
self-energy is computed by solving an auxiliary quantum 
impurity problem embedded in a self-consistent medium. 
We used the numerically exact continuous time quantum 
Monte Carlo method [12] . The on-site coulomb repul- 
sion on Fe-3d bands is estimated to be 4eV [5] and the 
Hunds coupling of J = 0.7 eV. The temperature is fixed 
at 116 K. 

The local spectral function A(ou) = A(k, uj) at tem- 
perature T = 116K is shown in Fig. [3^ together with the 
corresponding LDA density of states. The DMFT ap- 
proach predicts a renormalized low energy band with a 
fraction of the original width (Z ~ 0.2 — 0.3) while most 
of the weight is transferred into a broad Hubbard band at 
the binding energy ~ — 4eV. The system remains metal- 
lic at finite temperatures although very bad metal with 
the scattering rate at the Fermi level as high as 0.4 eV 
at 116 K. With decreasing temperature the Fermi sur- 
face is shrinking and the semiconducting gap is likely 
to open at zero temperature. Indeed slightly enhanced 
Coulomb repulsion (U = 4.5 eV) opens the gap even at 
room temperature. The correlation enhanced splitting 
between different orbitals leads to separation of bands 
into those that act as fully filled or fully empty bands 
at low energy. The high energy Hubbard bands are only 
weakly effected by such splitting. The Coulomb repulsion 
thus strongly reduces carrier density and pushes the par- 
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FIG. 3: a) DMFT density of states compared with LDA DOS. 
b) Orbitally resolved Fe-3d density of states within LDA and 
DMFT. 




FIG. 4: Momentum resolved spectral function A4d(k, uj) 
within DMFT (color coding) together with the LDA band 
structure (dashed lines). The upper panel corresponds to the 
undoped parent compound while the lower panel shows the 
10% electron doped compound in the virtual crystal approx- 
imation. 



ent compound on the verge of the transition between a 
bad metal and a bad semiconductor. At the same time, 
the localization of electrons leads to local moment for- 
mation and enhancement of the spin susceptibility in the 
doped compound. 

Many of the unconventional superconductors are 
known to have a very simple low energy band structure. 
For example, in the copper oxides, a single band is cross- 
ing the Fermi level. Similarly, the Fermi surface of the 
Na doped cobaltates [13 has primarily a\ Q single sheet 
Fermi surface. The situation is very different in LaOFeAs 
within LDA. As can be seen in Fig. ^Sjp LDA predicts 
that all five Fe 3d orbitals have finite weight close to or 
at the Fermi level. The situation is not simplified when 
the Coulomb correlation is accounted for. The spectral 
weight splits into high energy incoherent part and low 
energy part, which is very asymmetric and considerably 
reduced due to proximity to the semiconducting state. 
Upon doping, the quasiparticle peaks move to the Fermi 
level, the scattering rate is reduced, and the system be- 
comes better conductor. Experimentally, doping leads 
to the superconducting ground state at low temperature, 
which is likely to be of unconventional origin. The cooper 
pairs are likely to be formed out of composite singlets of 
spin and orbital degrees of freedom. 



Fig. [4^i shows momentum resolved spectral function 
A(k^co)LL in color coding together with the LDA 
bands in the energy range between — 5 eV and 2.5 eV and 
momentum dispersion in the high symmetry directions of 
the first Brillouin zone. The higher energy band structure 
at uj > 2eV and uj < — 5eV is not considerably different 
in the two approaches. In the intermediate frequency re- 
gion between — 1.5 eV and 1.5 eV a depletion of the spec- 
tral weight is apparent in DMFT approach. While LDA 
predicts a large number of bands in this region, DMFT 
redistributes most of this weight further away from the 
Fermi level to the range between — 2eV and — 4eV. Con- 
sequently a set of states with large scattering rate and 
short lifetime is predicted by DMFT. Finally the low en- 
ergy part of the spectra is considerably modified when 
Coulomb correlation is taken into account. The hole 
pockets around F point are barely identifiable and the 
electron pockets centered at M and A almost disappear 
- they are only touching the Fermi level and the weigh 
below the Fermi level is only due to large scattering rate. 
Thus the semiconducting type of gap is opening at the 
Fermi level. The correlations thus enhance the crystal 
field splittings which ultimately lead to a state being a 
band-semiconductor acting as bad metal at finite tern- 
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peratures. 

The situation is different in the doped compound (see 
Fig. [4)3) where the electron pockets clearly cross the 
Fermi level. The band velocity and effective mass are 
considerably enhanced (3 — 5-times) while the scattering 
rate still remains large. Finally, the hole pockets around 
T remain highly scattered. 
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FIG. 5: Optical conductivity of LaOFeAs within DMFT: 
in-plane and out-of plane. 



itinerant approach might be adequate. This will involve 
a generalization of the spin-fluctuation exchange mecha- 
nism to a multiorbital situation. In adition the pairing 
interaction might involve not just spin but also orbital 
fluctuations within the d subshell. 

In conclusion, we studied the band structure of newly 
discovered superconductor LaOi-^F^FeAs, and we pre- 
dict the orbital and momentum resolved spectral func- 
tion and optical conductivity of the compound. Density 
functional theory predicts that a set of Fe 3d bands are 
crossing the Fermi level with no clear splitting into e g 
and tig manifold. The coulomb correlations among the 
six electrons in the set of five Fe-3d orbitals is strong 
enough to push the compound close to the metal insula- 
tor transition. The correlation enhance the crystal field 
splitting among the d orbitals which in turn increase cor- 
relations. At temperature T = 116 K studied here, the 
parent compound LaOFeAs is still metallic although a 
bad metal with very large scattering rate and strongly re- 
duced number of carriers. Doping the parent compound 
leads to electron pockets centered at M and A point in 
momentum resolved spectral function. 



We also compute optical conductivity of LaOFeAs 
within the DMFT in both xy plane and along the z 
direction. Considerable anisotropy can be identified in 
Fig. [4Jd. The Drude peak, which is a hallmark of metal- 
licity and Fermi liquid state, is absent because of proxim- 
ity to the semiconducting state. Our calculations predict 
small number of carriers and large scattering rate, due 
to partial localization of carriers in this compound, and 
consequently low conductivity. One can identify a shoul- 
der in optical conductivity close to 2eV. It corresponds 
to transitions between Fe-3d and As 4p states. The large 
peak at 4eV is mostly due to transitions within the d 
shell from the lower Hubbard band to the quasiparticle 
peak above the Fermi level. 

In this work we used the single site Dynamical Mean 
Field Theory which approximates the full momentum de- 
pendent dynamical self-energy by its purely local compo- 
nent (see Eq. [T]) which is still a 10 x 10 matrix in 
orbital space of Fe-d orbitals. While this approximation 
successfully describes physics of many / and d electron 
systems, it does not capture the low energy physics of 
cuprate superconductors where momentum space differ- 
entiation leads to Fermi surface of arc-shape. The clus- 
ter extensions of DMFT can capture some of the reach 
physics of cuprates. We can not exclude the possibility 
that the low energy band structure of doped LaOFeAs 
will also require similar non-local correlations, and we 
leave this question open for future work. 

Finally we comment on possible routes to understand 
the superconductivity near a bad semiconductor in this 
compound. Since the Coulomb correlation is below the 
critical U for this system to undergo a Mott transition, an 
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